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New results of 300 hours of operation of the Tunka array are presented. An improved 
parametrization of the Cherenkov light lateral distribution function (LDF), based on 
CORSIKA Monte Carlo simulations and the experiment QUEST, has been used for the 
reconstruction of EAS parameters. The corrected energy spectrum in the knee region is 
obtained. The mean depth of the EAS maximum has been derived both from the analysis 
of LDF steepness and the FWHM of Cerenkov light pulse. The mean mass composition 
around the knee is estimated. 
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Experimental Differential Energy Spectrum. 

The TUNKA EAS Cerenkov array is located in Tunka Valley, at an altitude of 
675 m a.s.l., and was described in [1]. The new fitting function (LDF) for the EAS 
Cerenkov light lateral distribution, derived from CORSIKA code simulation [2], 
has been applied to TUNKA data. The primary energy Eq \TeV] has been obtained 
from the measured Cherenkov light flux at a distance 175 m from the shower core 
Qi75[photon ■ cm-'^ ■ eV''^] with CORSIKA sumulated relation: Eq = 370 • Q^yf 
The absolute energy calibration is based on the results obtained with the QUEST 
experiment [3]. A Monte Carlo simulation of the experiment has shown that the 
energy resolution is better than 18%. 

The spectrum is derived from data taken in 300 hours, spread over 51 clear 
moonless nights, with a trigger rate of about 1.8 Hz. To construct a spectrum with 
energy threshold 10^^ eV, showers with zenith angles 9 < 25° and a core position 
inside the geometrical area of the array have been selected. For the range from 
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Fig. 1. Differential energy spectrum. 

6- 10^^ to 10^^ eV, events with zenith angles 9 < 12° faUing inside a 5 times smaUer 
area around the array center have been selected. 

Estimation of Primary Mass Composition 

Lateral and time distributions of EAS Cherenkov light provide two independent 
methods to estimate the maximum depth. The simulation shows, that the LDF 
steepness P = (5(100)/Q(200) is related to the linear distance (in [km]) from the 
array to the EAS maximum position: H^nax = 10.62 - 0.12 ■ (P + 2.73)^. 

The Cherenkov light pulse FWHM [ns] at distances larger than 200 m from 
the EAS axis is related to the relative position of the EAS maximum by AX = 
Xo/cos9 — Xjnax [g/cHi^], whcrc Xq is the total depth of the atmosphere and 9 is 
the zenith angle of the shower. This relation depends only on the distance to the 
EAS axis. For example, for distance of 250 m: AX = 1677+1006 • logio{FWHM). 
This method gives a better theoretical accuracy, than the first one. Moreover, the 
Xmax estimate does not depend on assumptions about the primary nucleus. 

Figure 2 presents the mean depth of the EAS maximum, derived with the two 
methods described above, as a function of primary energy. It is seen that the thresh- 
old of the FWHM method is higher than that of the LDF steepness method, but 
the mean depths, obtained with the two different methods are in good agreement. 

The mean values of X^ax from fig. 2 can be easily transformed to the mean 
logarithmic mass < InA > of primary particles. Figure 3 represents the result. 
A slight correction, derived from MC simulations of the experiment (assuming a 
4-group mass composition, p:He:CNO:Fe — 0.3:0.3:0.2:0.2) has been applied. Ac- 
cording to these data the mass composition has almost no energy dependence in 
the range from 10^^ to 10^^ eV and is compatible with the hypothetical composi- 
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Fig. 2. Mean Xmax vs primary energy Eq. Fig. 3. Mean < InA > vs primary energy. 
Curves are CORSIKA/QGSJET simulations. 



tion with < InA >— 1.75 mentioned above. At energies above 10^^ eV, however, a 
steep increase of the average mass is observed. A similar increase in average mass 
is observed in almost all experiments. Points of the present work coincide with 
KASCADE data [4] at the energy of the knee. 
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